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1 Introduction 


Einstein, generalising the speeial theory of relativity to inelude gravita¬ 
tional foree gave the beautiful theory of general relativity wherein the spaee- 
time geometry took the eenter stage with its eurvature personifying gravi¬ 
tation. As a eovariant theory General relativity is eompletely independent 
of observers and their state of motion. Thus it is always interesting to 
look for adapting systems of eoordinates wherein one ean find some new 
interpretation that eould help in understanding the physieal behaviour of 
a system better. In faet it is well known that in the full 4 dim. spaee- 
time theory General Ralativity dispensed away with the Newtonian notion 
of inertial forees. However in 1988, Abramowiez Carter and Lasota showed 
that if one introduees the so ealled optieal geometry with a 3+1 splitting of 
spaee-time using a eonformal slieing, then one ean reintroduee the notion 
of Newtonian forees, whieh ean be useful in elarifying the origin of eer- 
tain GTR effeets that were known in the 4-dim. theory. For a free partiele 
only under the infiuenee of gravity the equations of motion are given by 
the geodesies of the given spaeetime manifold deseribing the gravitational 
field. Mueh as one admires the language of geometry the question is, does 
one get all the information inherent in the system within this framework 
or eould there be something overlooked or missed in the interpretation. 
Further, the language of forees has indeed been very useful in deseribing 
the other interaetions and thus it eould be useful to bring it baek into the 
realms of General Relativity, whieh might foeus the non-Newtonian fea¬ 
tures inherent in the geometry more explieitly. This approaeh would par- 
tieularly be of help in studying physieal aspeets very elose to ultra eompaet 
objeets and blaekholes, a region whieh is never available for the weak field 
Newtonian physies. 
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2 Formalism 


As the idea is to bring in Newtonian language into a geometrie theory of 
spaeetime, one needs to sliee the 4-spaee into a (3 spaee + time) strueture 
and look at the features on the absolute 3 spaee so obtained. In faet, sueh 
a 3+1 split of spaeetime is nothing new in general relativity, as, long ago 
Arnowitt, Deser and Misner (1962) introdueed sueh a seheme while look¬ 
ing for a method to give a Hamiltonian deseriptin of general relativity. As 
Misner, Throne and Wheeler (1972) mention, “The slicing of spacetime into 
a one parameter family of space-like hypersurfaces is called for, not only 
by the analysis of the dynamics along the way, but also by the boundary 
conditions as they pose themselves in any action principle of the form - give 
the 3-geometries on the two faces of a sandwich of spacetime and ac^ust 
the 4-geometry in between to extremize the action". 


In faet, sueh a proeedure of studying the dynamies effeetively paved the 
way for setting up numerieal methods to study evolution equations of the 
Cauehy data given on an initial hyper surfaee, and thus beeame a stan¬ 
dard proeedure for numerieal relativity (Seidel (1996) and York (1979)). 


Instead of a fully dynamieal system, suppose one has a stationary sys¬ 
tem wherein a time-like Killing veetor exists then one ean get a lower di¬ 
mensional quotient spaee through an isometry group aetion and one ean 
study eertain dynamieal features within a given geometrieal baekground. 
Abramowiez, Carter and Lasota (1988, hereafter referred to as ACL) used 
sueh a preseription with a eonformal resealing faetor and showed that 
one ean indeed obtain a 3+1 splitting wherein the 3-spaee is the quotient 
spaee obtained from the aetion of the time-like Killing veetor and the met- 
rie eonformal to the spatial geometry of the original four-spaee. As they 
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realised the most signifieant feature of a eonformal reslieing was that the 
normally geometrieal geodesie equation for a test partiele would separate 
into language of Newtonian forees wherein one ean direetly interpret terms 
as graviational, eentrifugal and Coriolis aeeelerations. 


Abramowiez, Nurowski and Wex (1993, hereafter referred to as ANW) later 
gave a eovariant approaeh to this formulation whieh does not depend upon 
any partieular symmetry and is as follows: 


In the given spaeetime manifold M with the metrie 

ds‘^ = gijdx'dx^ ( 1 ) 

introduee a eongruenee of world lines whieh is globally orthogonal to t = 
eonst., hypersurfaee whieh ensures that the vortieity of the eongruenee to 
be zero. In faet, Bardeen (1972) adopted sueh a eongruenee in axisym- 
metrie, stationary spaeetimes defining what are ealled loeally non-rotating 
observers or zero angular momenturm observers (ZAMO). The advantage 
of having sueh a eongruenee is that these loeal observers ‘rotate with the 
geometry’ and the eonneeting veetors between two sueh observers with ad- 
jaeent trajeetories do not preeess with respeet to Fermi-Walker transport. 
Denoting sueh a veetor field by rd {uird = — 1 time-like) it ean be verified 
that the eorresponding four-aeeeleration is proportional to the gradient of 
a sealar potential. 

n^^kni = Vicj ); = 0 ( 2 ) 

Though the veetor field rd is not uniquely determined by (2), loeally eaeh 
partieular ehoiee of rd uniquely defines a foliation of the spaeetime into 
sliees eaeh of whieh represents spaee at a partieular instant of time, whose 
geometry is given by 

hik = Qik + UiUk ; h\ = 5\ + rduk (3) 
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(2) also ensures that the speeial observers n* see no ehange in the poten¬ 
tial as their proper time passes by and thus have fixed positions that help 
them in distinguishing between different ‘inertial forees’. 


Consider a partiele of rest mass mo and four-veloeity U\ whieh ean be 
expressed as 

C* = 7 {rC' + rr*) (4) 

wherein r* is the unit tangent veetor (spaee-like) orthogonal to n* and par¬ 
allel to the 3-veloeity v of the partiele in the 3-spaee (Lorentz speed) and 
7 the Lorentz faetor ) • The four-aeeeleration of the partiele a* may 

now be obtained through direet eomputation (Abramowiez, 1993) 

Ofc = = - 7 ^Vfc 0 -I- 7^r -t iUk) 

9 9 

+72^VVjTfc -t (r7)-rfc -t 7nfc 

Let us eonsider the motion of the partiele in a eireular orbit in a general 
stationary axisymmetrie spaeetime. From the given symmetries, there ex¬ 
ist two Killing veetors 7* the time-like having open trajeetories, and Q the 
spaee-like with elosed trajeetories. If the partiele has a eonstant angu¬ 
lar veloeity Vt as measured by the stationary observer at infinity then its 
four-veloeity may be expressed as 

u^ = A (rf + flC*) (6) 


A being the redshift faetor 


42 


<r]r] > < 77C > < CC > 


(7) 


In terms of the Killing veetors one ean express n* and 0 eonsistently as 


n* = (7* + uC) , uj = - <r],C > / < C,C 


> 


( 8 ) 


and 


= —in 
2 


<T],T]> -2uj <C,V> < C, C > 


(9) 


5 



From (4) and (6) one can evalutate the particle speed V to be: 

Vt^ = e'^(Q-aj)C ( 10 ) 

Using now the ACL approach of conformal rescaling of the 3-metric, one 
can define the projected metric and the vectors 

hij = hij , f* = e'^r* (11) 

such that the acceleration may now be written as 

ak ■■= -Vfc0 + 7V (n^ViTk + T^ViUk) + {-fVfrVifk (12) 

as the last two terms in (5) become zero for a particle with constant speed, 
and conserved energy. V in (12) refers to the covariant derivative with re¬ 
spect to the metric hij. As may be seen, the acceleration is made up of 
three distinct terms, (i) gradient of a scalar potential, (ii) a term propor¬ 
tional to V, and (iii) one proportional to The first and the third terms 
may be immediately recognised as the gravitational and centrifugal accel¬ 
erations. Further, as n* and r* are parallel to the Killing vector 7* and C 
respectively from the equation for Lie derivative, one has 

LnT = rhWiTk + TiVku" = 0 (13) 


Using this alongwith the fact that n* and r* are orthogonal, the second term 
of (12) may be written as 
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rh (ViTfc - VfcT*) 


(14) 


which represents the Lense-Thirring effect of the inertial drag and thus 
identified as the generalisation of the Coriolis acceleration. 


If the general axisymmetric and stationary spacetime is represented by the 
metric 

ds'^ = gttdf + 2gtj>dtd(t) + gj>j,d(f + grdr^ + geedO"^ ( 15 ) 
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with gijS being functins of r and 9 only then the total foree aeting on a par- 
tiele in eirular orbit with eonstant speed Q may be expressed as (Prasanna 
(1997)) 


wherein 


and 


with 


and 


F, := (Gr), + (Go), + (G/), 


(Gr), : 
(Go), : 


-Vi(j) = ^di [in - gtt9<p,p) /g 4 ><p \} 


= (VfcTj - VjTfc) 


= —a 




(G/)_ = (7F)2f^Vfcf, 

= -A^{Q - [in [gl^\ [g^^ - gug^^)] } 


$ = 


gtt ^ g4'4> ’ ^ 9tq/ 94'4> 


A^ = - 


Qtt + 9<i><t> 


-1 


(16) 


(17) 


(18) 


(19) 


In order to understand this splitting of the total foree, let us eonsider the 
simple ease of statie spaeetimes (Abramowiez and Prasanna (1990)). Then, 
by definition the Coriolis term would be absent, and the foree aeting on a 
test partiele of mass mo and 3-momentum p" is given by (ACL) 

2 

~ ~ Tfl ~ 

mofa = ( 20 ) 


wherein = moCyf^ and thus eonsistent with what one gets from (15) in 
the projeeted 3-spaee. It is now elear that in the statie spaee time, the 
trajeetories of photons (rest mass zero partiele) are given by the eurves 


= 0 


( 21 ) 


whieh by definition are geodesies of this 3-spaee. Thus, one finds that the 
‘null trajetories’ of the 4-spaee projeet onto the geodesies of the quotient 
spaee obtained through eonformal slieing and indieates that partieles on 
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these trajeetories do not experienee any ‘foree’ in the Newtonian sense. 


It is for this reason that ACL ealled this slieing as optieal referenee geome¬ 
try meaning that the null lines projeet onto straight lines in the Euelidean 
sense. However, as was seen later this is true only for statie spaeetimes, 
as rotation would influenee the photons differently for prograde and ret¬ 
rograde motion. In statie spaeetimes, one finds that partieles aeted on 
by forees other than gravity would deviate from the geodesies of the quo¬ 
tient 3-spaee obtained by eonformal slieing, as these geodesies behave like 
straight lines of Newtonian geometry and thus follow Newtonian laws of 
motion. 


The most important aspeet of the frmalism is that while the partiele kine- 
maties is expressed in Newtonian language, the general relativistie effeets 
are all present and thus help in understanding results whieh were earlier 
known in general relativistie analysis but were hidden in the geometry and 
not aeeessible for visualisation. 


3 Specific Applications 

We start from the simplest applieation of the methodology outlined above to 
study the partiele kinematies in the statie spaeetimes, taking the Sehwarzsehild 
geometry as the first example (Abramowiez and Prasanna, 1990; hereafter 
referred to as AP). 


The metrie as expressed in the usual eoordinates, 

-I- {dO^ + sin^ Odcp^'^ (22) 
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would yield for the gravitational and eentrifugal aeeelerations, aeting on a 
test partiele in eiruelar orbit, the expressions: 




-1 


and 




(23) 


(24) 


From (24) it is elear that while {Gr) and (C/) are in opposite direetions upto 
r = 3m from infinity for r < 3m they both have the same direetion. As gravi¬ 
tational aeeeleration is always undireetional, it is elear that the eentrifugal 
aeeeleration reverses its sign at r = 3m. In faet, if one looks at the photon 
effeetive potential in the Sehwarzsehild geometry, one finds that r = 3m is 
the loeation of the maximum and thus eorresponds to unstable eireular 
orbit. What has been notieed now is that this null line is the straight line 
path for the photon in the quotient spaee and thus eorresponds to a loea¬ 
tion at whieh the eentrifugal aeeeleration is zero, and on either sides the 
eentrifugal foree aets in opposite direetions. 


As shown in AP, this feature has important kinematieal implieations in the 

study of aeeretion fiows near ultra eompaet objeets (blaek holes) like the 

Rayleigh eeiterion for stability of fiow turns out to be 

2m(r^ —2mr^) 3m\ ^ ^ 

(r3 - i‘^r -t \ r J dr 

This means ^ > 0 for r > 3m and < 0 for r < 3m, indieating that for r < 3m 
the angular momentum has to be adveeted inwards for stability. In faet, 
this result elearly explained the findings of Anderson and Lemos (1988), 
who had obtained inward adveetion of angular momentum very elose to 
blaek holes. 


Another important implieation of eentrifugal reversal is borne out in the 
evaluation of elliptieity of slowly rotating fiuid eonfiguration represented 
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by a sequence of quasi-stationary solutions, with decreasing radii, keeping 
the mass and angular momentum conserved. Chandrasekhar and Miller 
(1974) had considered this problem and found that the ellipticity instead 
of increasing continuosly, reached a maximum, a result which they had 
attributed to frame dragging of rotatintg systems. However, after the dis¬ 
covery of centrifugal reversal, Abramowicz and Miller (1990) analysed the 
equilibrium configuration in general relativity using the ‘inertial forces’ ap¬ 
proach and found that the equilibrium demands 

02 _ GM 

“fc - R3 

which to the lowest order in Q gave the ellipticity function to be 


e{R) 


125 _ 3 ■ 

I ~ 



(27) 


exhibiting a maximum at i? = 3. 


Though they had obtained the ellipticity maximum, the lacuna in their 
approach was that they had used only the Schwarzschild exterior geometry 
for evaluating the inertial forces. Prasanna and coworkers (A. Gupta, S. 
Iyer and A.R. Prasanna (1996)) reconsidered this problem starting from the 
general conservation laws expressed in the 3+1 formalism and using the 
Hartle-Thorne approximation solution for the interior of a slowly rotating 
body. Defining the ellipticity through force balance equations at the pole 
and equator 


e2 = l-^ 


gej 


e ^ — 


(28) 
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[Fge is the gravitational force at the equator and Fgp at the poles) they found 
the maximum to occur at r = 5.4m, which is much closer to Chandrasekhar- 
Miller result of maximum occurring at r = 5m. 
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Prasanna (1991) had considered another statie spaeetime viz. Ernst spaee- 
time whieh represents the external field of a blaekhole immersed in a uni¬ 
form magnetie field. The photon effeetive potential in this metrie has two 
extrema, the maxima eorresponding to the unstable orbit very elose to 
r = 3m, while a minima loeated far away, governed by the strength of 
the magnetie field and yielding a stable photon orbit. As eentrifugal foree 
would go to zero at both these loeations, it appears that the Newtonian 
orbits are possible only within the region bounded by these two points of 
extrema. 


We next eonsider the situation in the Kerr geometry whieh is aetually sup¬ 
posed to be the spaeetime exterior to a rotating blaekhole represented by 


the metrie 


ds^ = {l - dt^ - ^ smHdtd(j) + 


with 


■ sin^ 6d(j)‘^ 

A = — 2mr + 

S = r‘^ + a? cos^ 9 


(29) 


B = 


(r^ -I- a^) — Aa^ sin^ 9 


Using the formluae presented in (16)-(18), one ean get (Prasanna (1997)) 

m (a^A + + a^r‘^ — 2ma?r) 


{GR)r = 


rA (r^ -|- a^r -|- 2ma^) 


(30) 


{Cf)r 


— (U — ujY [r^ 


r^A 


1 


— 3mr^ + (r^ — 3mr^ -|- 6m^r — 2ma?‘) 

- 02 (r2 + a?) - (1 - naf 


(31) 


{C0)r = 


2ma{Q — u) (3r^ -|- a^) 


r (r^ + a?r + 2ma?) 1 — (r^ -|- a?) — (1 — Ua)" 


(32) 


One ean immediately see the differenee in the nature of eentrifugal and 
Coriolis forees, whereas the Coriolis depends on the eoupling of the angu¬ 
lar momentum of the eentral souree with that of the partiele a{Vt — u), the 
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centrifugal can go to zero at different locations solely depending upon ‘a 
due to the zeros of the quintic equation 

r® — 3mr^ + (r^ — 3mr‘^ + 6m^r — 2ma^'^ = 0 (33) 



0.2 0.4 0.6 0.8 1 


Fig. 1. Location of the point where Cfg = 0, (-), retrograde photon 

orbit (-) and of the prograde photon orbit (-■ -■ -■) for different values of 

a. 


It may be easily seen that for 0 < a < 1 the equation (32) at best can have 
only three real roots of which, one is always definitely outside the event 
horizon (Iyer and Prasanna (1993)) as depicted in Fig. 1. However, as also 
shown in this figure this location does not coincide with the location of the 
unstable photon orbit, prograde or retrograde. The centrifugal force van¬ 
ishes at a location between the two photon orbits and for the case a = 0, 
they all coincide at r = 3m. Unfortunately, the direct link between the 
unstable photon orbit and the centrifugal force reversal, depicted in static 
spacetimes do not find a parallel in stationary spacetime. Rotation does 
indeed bring in some new features of which the frame dragging is the most 
important one. 
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Figs. 2 and 3 show the nature of the eentrifugal and Coriolis forees at 
the loeation of retrograde and direet photon orbits as a funetion of Q for 
a = 0.5. The first impression that one gets is that these forees ehange sign 
for different values of Q aeross the asymptotes. However, one has to eheek 
that the asymptotes are eaused by the infinity of the redshift faetor at 
the roots of the equation 

+ (?« = 0 (34) 

= cu ± \/u‘^ - 9tt/9<t><t> (35) 

Henee the only portion of the plots whieh is meaningful is the region eorre- 
sponding to the values of H; H < H+. As may be seen in this region the 
eentrifugal is positive along the retrograde photon orbit [rph+) ane negative 
along the direet photon orbit [rph—) as it should be aeeording to Fig. 1. On 
the other hand, the eoriolis ehanges sign in both eases at the point H = c.;, 
wherein eentrifugal also vanishes beeause of the faet that the angular ve- 
loeity is just equal to the dragging of inertial frames, by the spaeetime due 
to the rotation of the eentral objeet. 


On the other hand, we have seen that there are loeations in the given 

spaeetime, wherein the eentrifugal foree is zero but Q ^ u and thus the 

Coriolis is non-zero. In view of this, Prasanna (1997) defined an index of 

referenee ealled the ‘Cumulative Drag Index’ as defined by the ratio 

{Co-_^ 

{Co + Gr) 

whieh eould eharaeterise purely the rotational feature of a spaeetime through 
its infiuenee on a partiele in eireular orbit at the loeation where the een¬ 
trifugal foree is zero. 
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Fig. 2. Centrifugal (-) and Coriolis (-) force plots for a = 0.5, along 

the retrograde photon orbit for — 1 < < 1. 



Fig. 3. Same as (2) at the prograde photon orbit. 









Fig. 4 shows the plot of C as a function of fl, for a fixed a and R. As 
may be seen, there are two zeros and two infinities for the function. As 
a > 0, > 0 represents the co-rotating particles and < 0 represents the 

counter-rotating particles. 



Fig. 4. Cumulative drag index C for a = 0.5, R = 2.9445 as a function of if 

(-1 < n < 1 ). 

As the orbit chosen is the one where the centrifugal force is identically 
zero the infinities of C refer to the trajectories along which the total force 
acting on the particle is zero as given by the definition of C. Fig. 5 shows 
the index C for three distinct cases of a = 0.1, 0.5 and 1 at the respective 
locations where Cf = 0. 


It is interesting to note that as a increases, the co-rotating particles have to 
decrease their angular velocity Q very little to stay in equilibrium, whereas 
the counter-rotating ones have to increase their Q much more, as depicted 
explicitly in Table 1. 


15 



Cdl 



Fig. 5. C for three different values of a showing differenee in values for 
equilibrium orbit {Co = —Gr). 


Table 1 


a 

R 




u 

0.1 

2.9978 

0.189022 

-0.1964 

0.0074 

0.0074 

0.5 

2.9445 

0.180094 

-0.21965 

0.0395 

0.0374 

1.0 

2.7830 

0.17722 

-0.27452 

0.097 

0.076 


However, it may be seen that, the exeess adjustment on the part of eounter- 
rotating partieles is essentially due to the dragging of inertial frames u 
whieh is always in the direetion of rotation of the eentral souree. Table 
1 gives the values of u for given a and R, whieh matehes quite elosely 
with the differenee in the eorresponding Q for given a and R, between eo- 
and eounter-rotating partieles. Thus one finds that using the language of 
forees within the framework of general relatvity, one ean exaetly quantify 
the spaeetime eurvature effeets whieh otherwise is hidden in the geometry. 


There have been other approaehes to deseribe partiele motion in general 
relativity a general survey of whieh was provided by Bini et al. (1997). 
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However, the eonformal slieing as deseribed above helps direetly in un¬ 
derstanding eertain physieal features of fluid flow eonflgurations near ex- 
tremaly eompaet objeets and also visualise the inertial drag more elearly. 
It would indeed be very niee if one ean use this in the eontext of ADM 
formalism for Hamiltonian dynamies. This eould be of direet eonsequenee 
to understand the dynamieal eollapse of fluid eonflgurations, through nu- 
merieal approaeh. It is indeed very neeessary to follow the eollapse, par- 
tieularly to see whether the eentrifugal reversal that is inherent would 
produee the strange behaviour of the elliptieity as evideneed in the work 
of Chandrasekhar and Miller. A sudden ehange in elliptieity eould indeed 
produee notieeable signature on the emission of gravitational radiation 
from a eollapsing star in the final stages before it gets into a blaekhole. 
This is an open problem whieh eould be of great interest for elassieal grav¬ 
itation studies. 
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